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Abstract 



The energy dependence of the total kaon production cross sections in baryon baryon 
(N and A) collisions are studied in the resonance model, which is a relativistic, tree-level 
treatment. This study is the first attempt to complete a systematic, consistent investi- 
gation of the elementary kaon production reactions for both the pion baryon and baryon 
baryon reactions. Our model suggests that the magnitudes of the isospin-averaged total 
cross sections for the NN — ► NYK and AN — ► NYK (Y = A or E) reactions are almost 
equal at energies up to about 200 MeV above threshold. Bowever, the magnitudes for the 
AN reactions become about 6 times larger than those for the AW reactions at energies 
about 1 GeV above threshold. Furthermore, the magnitudes of the isospin-averaged total 
cross sections for the NN — ► AYK reactions turn out to be comparable to those for the 
JVJV — * NYK reactions at NN invariant collision energies about 3.1 GeV, and about 5 to 
10 times larger at NN invariant collision energies about 3.5 GeV. The microscopic cross 
sections are parametrized in all isospin channels necessary for the transport model studies 
of kaon production in heavy ion collisions. These cross sections are then applied in the 
relativistic transport model to study the sensitivity to the underlying elementary kaon 
production cross sections. 
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1 Introduction 



Particle production in heavy ion collisions is a unique tool for studying the properties of matter 
under extreme conditions such as high temperature and/or density [0-0]. In particular, pos- 
itive kaon, K + , production in heavy ion collisions at intermediate energies is one of the most 
promising methods to probe matter formed under such conditions in the central zone of the 
collisions [|[ §]. 

Because 7^ + -mesons have a long mean free path due to the small cross sections for scat- 
tering on nucleons, they can provide us almost undisturbed information about the central zone 
of heavy ion collisions, in which most of the i^ + -mesons are considered to be produced in sub- 
threshold heavy ion reactions. According to the theoretical investigations in Refs. [f7|- 1|1 1|| , the 
kaon yield is also very sensitive to the nuclear equation of state. Furthermore, kaon production 
in heavy ion collisions may provide clues concerning quark gluon plasma formation from the 
K/tt production ratio (JT21 |T^], chiral symmetry restoration [|TJ], [T5|, IB|, and the existence of 
density isomers [17| . Thus, many studies of kaon production in heavy ion collisions have been 
made experimentally and theoretically [|I1|-[|62|. 



The emphasis and focus of earlier theoretical studies was on the momentum- dependent 

m 



nuclear interactions j2C| [2 
effects of chiral symmetry restoration |TJ 
production cross sections 



22], effects of the rescattering and kaon mass renormalization 



and the medium modification of elementary kaon 
25). However, most of these theoretical investigations were per- 
fomed using the energy dependence of the total kaon production cross sections parametrized in 
a simple phenomenological way P, I 
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As a consequence, all these investigations involve a 
certain amount of ambiguity arising from such parametrizations. Furthermore, although recent 
investigations of strangeness production emphasize the importance of elementary processes such 
as the BiB 2 
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no investigations for these reactions 



B3YK and txB — ► YK reactions 
have been perfomed in a consistent manner even in a tree-level treatment. (Here -81,2,3 an d B 
stand for either the nucleon or the A, while Y stands for either the A or the S hyperon.) 

Although many issues need to be considered in order to achieve a better understanding 
of kaon production in heavy ion collisions, the one emphasized in this article is the elementary 
kaon production cross sections at the hadronic level as was studied in Refs. [ER 
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6J . Since most of the microscopic transport models assume that the relevant reactions 
are described by on-shell, binary interactions at the hadronic level, it should be straightforward 
to improve the parametrizations for the elementary kaon production cross sections, which are 
important inputs of the models. 

Recently, the COSY- 11 Collaboration |3^] measured the total cross section for the pp — > 
pAK + reaction at an energy 2 MeV above the reaction threshold. The result distinctly differs 
from those of the commonly used phenomenological parametrizations || pjfl . However, in 
practice, this discrepancy probably does not significantly influence the past results calculated 
using the parametrizations made without this new data point, because the magnitudes of the 
cross sections at this energy are very small for both the experiment and parametrizations 0, 



On the other hand, theoretical reanalyses in Refs. [^, [36], |37| using different parametriza- 
tions [3|[] for the iViV — > NY K cross sections, indicate that the secondary reactions, 7riV — > YK, 
give rather important contributions - comparable with those of the baryon baryon (N and A) 
induced reactions. This conclusion contradicts the previously accepted scenario, in which the 
dominant contributions for K + production in heavy ion collisions come from AiV and AA 
collisions [10. |39, HO, O, 42|. Thus, it is necessary to investigate further the elementary kaon 
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production cross sections for both the ttB — > YK and BiB 2 — > B 3 YK reactions in a consistent 
manner. Furthermore, all the conclusions drawn about the elementary reactions involving a 
A were based on the results obtained using the elementary cross sections for the nucleon, not 
for the A, except for the isospin difference. Thus, an explicit calculation for the cross sections 
involving a A should be performed, treating properly not only isospin degrees of freedom but 
also spin as was initially attempted in Ref. [42 using the results of the resonance model for the 
7rB — > YK reactions [[46] . 

For this purpose, we complete our systematic investigations for the elementary kaon 
production processes in all baryon baryon reactions, extending the studies made so far for the 
pion baryon and proton proton reactions using the resonance model |44 , [45| . The treatment of 
the present study is a tree-level, using empirical branching ratios for the relevant resonances. It 
is relativistic, and incorporates the pion baryon, and the baryon baryon reactions in a consistent 
manner. Thus, most of the cut-off parameters, coupling constants and form factors have been 
fixed already in pion baryon and proton proton reactions p4l f46|. Furthermore, we parametrize 



the energy dependence of these total kaon production cross sections in baryon baryon collisions 
for all the isospin channels needed in the simulation codes. We note that the total cross 
sections for the pp — > pAK + reaction recently measured at energies about 50 and 150 MeV 
above threshold [4? show an excellent agreement with the predicted results of the resonance 
model 
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However, because the treatment in the model does not include the final state 
interactions, our results cannot be expected to describe well the near-threshold behaviour. A 
more rigorous treatment of the resonances consistent with the unitary condition was studied 
by Fuester and Mosel JKJ, but such an approach does not seem to be suitable for the present 
purpose, and beyond the scope of the present study. 

The organization of this article is as follows. In Sect. 2 we will explain the resonance model 
in detail. Effective Lagrangian densities at the hadronic level together with the experimental 
data for the model will be described. Numerical results for the total cross sections will be 
presented in Sec. 3. In Sect. 4 we give all the parametrizations for the energy dependence of 
the total cross sections calculated in the model. In Sec. 5 these cross sections are used in 
the relativistic transport model to study the sensitivity of the transport model predictions to 
the underlying elementary kaon production cross sections. Discussions and summary will be 
given in Sect. 6. Finally, in the Appendix, the relations between the branching ratios of the 
resonances and the coupling constants relevant for the model will be given. In addition, we will 
supplement the cross section relations for the ttN — > AK reaction which were not mentioned 
in Ref. EH . 



2 Resonance model 

In Fig. |l|, we show the Feynman diagrams relevant for kaon production in the resonance model. 
R stands for the baryon resonances which are responsible for the kaon and hyperon pair pro- 
duction, with their masses up to about 2 GeV. Assumptions and approximations of the model 
are @,g§,g6J: 



1. Resonances which are experimentally observed to decay to hyperon and kaon are included 
in the model. A kaon is always produced from these resonances simultaneously with a 
hyperon. In the study of the ttB — ► YK reactions, we have also investigated the non- 
resonant contributions, including an effective t-channel i^*-meson exchange. It turned 
out that this non-resonant contribution was small in order to reproduce the experimental 
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data, once the relevant s-channel resonances were included in the model [46]. Thus, 
kaon exchange which is included in Refs. ER BO, pTl |38], O and would introduce extra 



new parameters in the present treatment, is not included in the present calculation. 
Instead, other meson exchanges, namely rj and p meson exchanges, are included, and 
they compensate the contribution of kaon exchange introduced in other models. 

2. Mesons exchanged are restricted to those observed in the decay channels of the adopted 
resonances. (See also Fig. [I].) We note that there is a discussion in Ref. of whether 
there is a possibility to settle what kind of meson exchange is responsible for the pp — > 
pAK + reaction. 

3. Those processes in which the exchanged pion can be on-shell are excluded, because on- 
shell pions are usually included as secondary processes in microscopic transport models 
through processes such as ttN <-> A and ttB — > YK. Thus, the AA — > NYK reactions, 
which are possible only through 7r-exchange in the model, are not included. Note that 
because of the assumption made in number 2, neither p-meson nor 77-meson exchanges 
will give contributions to this channel in the present treatment. 

4. Resonances are treated as elementary particles, except that their observed widths enter 
in the propagators. Because of this, p-exchange attached to the iV(1650) resonance is 
excluded kinematically, although this resonance has a branch in the N p channel |)3 . 



5. Any final state interactions which are usually important for describing the near-threshold 
behaviour are not included. Thus, our results should not be taken seriously near the 
reaction threshold, which is different from the present main purpose. 

It is worth noting that all values of the coupling constants squared relevant for the meson- 
baryon-(baryon resonance) vertices, Qmbr iM = KjVjP)-) can be determined from the exper- 
imental decay ratios, once the interaction Lagrangian densities and relevant form factors are 
specified. In Table p] we summarize the data for those resonances which are included in the 
model, and necessary to calculate the coupling constants. 

Effective Lagrangian densities relevant for evaluating the Feynman diagrams depicted in 
Fig. |] are: 



£>irNN 
£>ttNA 

£>ttNN(1650) 
£ttNN(1710) 
AtAW(1720) 

^ttVA(1920) 



-WttNnN^tN ■ 7T, 

_9 I NA ( ^ N _ g -, + ft^A» _ Q 



7T 



-g*NN(i6B0) (JV(1650)tW • 7? + iVf JV(1650) • vr 
-ig*NN(mo) (iV(1710)75riV ■ 7? + N l5 rN(1710) ■ rr ) , 
^ JVjV(1720) (A>(1720)fiV • dj + iVf^(1720) • dj 
ftrJVA(i920) ( AM(i920) ?N ■ + ^2^(1920) • d i 



7T 



7T. 



c 



ttAJV(1650) 



-^ttAaA m 7 5 /CA a 

.^AAr(i65o) ( ^^650)^^^ • dfi + A^tNilGbO) ■ 9 m tt 



C 



ttAJV(1710) 



g?rAJV(1710) ( N(1710) 2* f A" ■ + AfiV(1710) • d. 



7T 



(1) 

(2) 

(3) 
(4) 
(5) 

(6) 

(7) 
(8) 

(9) 
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Resonances included in the model. Confidence levels of the resonances are, iV(1650) * 

PI. Note that the A(1920) resonance is 



Table 1: 

* * *, jV(1710) * **, iV(1720) * * * * and A(1920) * ** 



treated as an effective resonance which represents all contributions of six resonances, A(1900), 
A(1905), A(1910), A(1920), A(1930) and A(1940). See Ref. @ for this effective treatment of 
the A(1920). 



Resonance (J F ) 


Width (MeV) 


Decay channel 


Branching ratio 


Adopted value 


N(1650)(±~) 


150 


Nn 


0.60 - 0.80 


0.700 






Nr] 


0.03 - 0.10 


0.065 






Att 


0.03 - 0.07 


0.050 






AK 


0.03 - 0.11 


0.070 


N(1710)(± + ) 


100 


Ntc 


0.10 - 0.20 


0.150 






Nr] 


0.20 - 0.40 


0.300 






Np 


0.05 - 0.25 


0.150 






Att 


0.10 - 0.25 


0.175 






AK 


0.05 - 0.25 


0.150 






Y.K 


0.02 - 0.10 


0.060 


N(1720) (| + ) 


150 


Nir 


0.10 - 0.20 


0.150 






Nr] 


0.02 - 0.06 


0.040 






Np 


0.70 - 0.85 


0.775 






Att 


0.05 - 0.15 


0.100 






AK 


0.03 - 0.10 


0.065 






Y.K 


0.02 - 0.05 


0.035 


A(1920)(| + ) 


200 


Nn 


0.05 - 0.20 


0.125 






Y,K 


0.01 - 0.03 


0.020 
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^Aivcmo) = -i^Aivcmo) ^(1720)75^ A" • t? + A" 75 lV(1720) - tt J , (10) 

C v nn = -ig V NNN~f 5 Nr), (11) 

Cr,AA = -^r,AAA M 7 5 A M ?7, (12) 

£„w(i650) = -9 V NN(i65Q) (n (1650) Nr/ + NN(lQ50)r] ) , (13) 

C V NN(i7io) = -^Arjv(mo) (iV(1710)75^77 + 7V75^(1710)?7 ) , (14) 

£^(1720) = ^ JVjV(1720) (^(1720)^ + ^(1720)9^) , (15) 

£ p iwv = -^iv (WtW • p M + -JL-Na^rN ■ , (16) 

£ p aa = ^aaA^A" • p„, (17) 

£ p na = -i^U^l^N + N^^AA-id^-d^), (18) 

£ P iwv(mo) = -^ATiv ( i7io) (iV(1710) 7 M riV • p M + iVyW(1710) • p M ) , (19) 

£^(1720) = -^ P JViv(i720) (iV"(1720) 75 7W • p M + #75^(1720) • p M ) , (20) 

Cranio) = -Qkan^o) (iV(1650)AK + XA7V(1650)) , (21) 

£*AJV(mc» = -^xAiv(mo) (iV(1710) 75 AX + XA 75 iV(1710)) , (22) 

/W(i720) = 9KAN{1720) (N^1720)Ad li K + (c^)AiV(1720)) , (23) 

£*Eiv(i7io) = -^xsiv(i7io) (iV(1710)7 5 r • + K ? • f 75 7V(1710)) , (24) 

£xsiv(i72o) = gJ " SJV(m0) f^(1720)f • ^fyif + (^X)S • f 7^(1720)) , (25) 

^EA(i920) = ^ SA(1920) ( A^(1920)^ ■ ifc^if + (d^Ryr, ■ 2^(1920)) . (26) 

In the above, the operators and A? are defined by 

^m, = E (i4^l? M )^' ( 27 ) 

€=±1,0 Z Z 
€=±1,0 Z 1 



with M, /x and M' being the third components of the isospin projections, and f the Pauli 
matrices. N, A^(1710), AT(1720) and A(1920) stand for the fields of the nucleon and the cor- 
responding baryon resonances. They are expressed by N = (p, n) , similarly for the nucleon 
resonances, and A(1920) = (A(1920)++, A(1920)+, A(1920)°, A(1920)~) in isospin space. The 
physical representations of the kaon field are, K T = (K + , K°) and K = (K~ , K°^j, respectively, 
where the superscript T means the transposition operation. They are defined as annihilating 
(creating) the physical particle (anti-particle) states. For the propagators iSf{p) of the spin 
1/2 and iG^ v (p) of the spin 3/2 resonances we use: 

iS F {p) = i— , , 29 

p z — m A + imi /M " 



5 



^(ri = V-m^r*» - (30) 



with 

ST ~ |w ~ ^(7Y - tV) ~ -^ff 



P^(p) = _( 7 . p + m ) 



(31) 



where m and stand for the mass and full decay width of the corresponding resonances. For 
the form factors, Fm{<1) (q is the momentum of meson, M), appearing in the meson-baryon- 
(baryon resonance) vertices, we use 

FlM = (aETf)" (32) 

where n = 1 for the 7r and 77- mesons and n = 2 for the p- meson, respectively, with Am being 
the cut-off parameter. Most of the form factors, coupling constants and cut-off parameters for 
the relevant meson-baryon-baryon and meson-baryon-(baryon resonance) vertices are adopted 
from Refs. [[0], (46| . In the Appendix, we give the relations between the branching ratios and 
the corresponding coupling constants squared, g MBR , which were calculated using the relevant 
Lagrangian densities. In the calculation, the form factors of Eq. fl32"l) are multiplied by the 
corresponding coupling constants, gMBR- For the coupling constants, <?a/aa (M = ir,r],p), 
which appeared for the first time in the present study, we use an SU(6) quark model result 
with the definition Eq. (|28|), <7maa = Sqmnn- In addition, we use the same value for the 
corresponding cut-off parameter as that for the nucleon. For the value of the cut-off parameter 
at the pNA vertex, we use 544 = 1300 x (920/2200) MeV, which is scaled the same amount as 
was necessary for the pNN vertex of the Bonn potential model (the Model I in Table B.l ||64|| ). 
The other quantities, coupling constants, cut-off parameters and form factors are, as far as 
possible, taken from the same version of the model []64 |. In addition, we use a value, k = 



fpNN 1 9pNN = 6.1, for the tensor coupling constant at the pNN vertex. We summarize in 
Table all values for the coupling constants and cut-off parameters used in the study. 



3 Numerical results 

In this study, we neglected all interference terms between the amplitudes. Thus, the relations 
given in this section for the cross sections are not always valid when the interference terms are 
included rigorously. 

3.1 NN -> NYK 

Recently, the total cross section for the pp — > pAK + reaction was measured by the COSY 11 



Collaboration at energies of a few MeV above the reaction threshold, and it gave a new 
constraint on the theoretical calculations and phenomenological parametrizations. However, 
we do not refit the parameters and coupling constants to this new data point for the following 
reasons. First, our treatment does not include the final state interactions nor interference terms, 
which will be important at energies very near threshold. Second, the total cross sections at 
energies up to about 10 MeV above the threshold are very small in magnitude for both the 
parametrization and the new data point. Thus, this new data point will not influence the 
calculation of kaon yield in heavy ion collisions (for instance see Fig. 9 in Ref . |}5| ) . 
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Table 2: Coupling constants and cut-off parameters used in the present study, k = 
fpNN 1 9pN n = 6.1 is used for the pNN tensor coupling. Note that the coupling constants 
relevant for A(1920), g n N A(i92o) an d 5 i ksa(i92o); are scaled multiplying by a factor 1.86 accord- 
ing to the effective treatment [pE6]| . See also the caption of Table [1]. 



vertex 


g 2 /^ 


cut-off (MeV) 




vertex 


g 2 /^ 


cut-off (MeV) 


nNN 


14.4 


1050 




ttAA 


3 2 x 14.4 


1050 


vriViV(1650) 


1.12 x 1CT 1 


800 




vriViV(1710) 


2.05 x 10- 1 


800 


ttAW(1720) 


4.13 x 1(T 3 


800 




?riVA(1920) 


1.13 x 10" 1 


500 


vrAiV(1650) 


7.19 x 10~ 2 


800 




vrAiV(1710) 


2.23 x lO" 3 


800 


vrAiV(1720) 


1.39 


800 




r]NN 


5.00 


2000 


r]AA 


3 2 x 5.00 


2000 




r]NN(165Q) 


3.37 x 10~ 2 


800 


rjNN(1710) 


2.31 


800 




r]NN(172Q) 


1.03 x 10" 1 


800 


pNN 


0.74 


920 




pNA 


19.0 


544 


pAA 


3 2 x 0.74 


920 




pN N(1710) 


3.61 x 10 +1 


800 


pNN (1720) 


1.43 x 10+ 2 


800 




KAN(1650) 


5.10 x 10~ 2 


800 


KAN (1710) 


3.78 


800 




KAN(1720) 


3.12 x 10" 1 


800 


KEN(1710) 


4.66 


800 




KEN(1720) 


2.99 x 10- 1 


800 


fQ]A(1920) 


3.08 x 10- 1 


500 











Here, we write down explicitly the amplitude and cross section formula for the pp —>■ pAK + 
reaction, as an example. The total amplitude for this reaction is given by 

M = A^(p(1650),7T°) + M(p(1710),tt°) + M(p(1720),7i°) 
+ M(p(lQ50),r]) + M(p(1710),rj) + M(p(1720) , rj) 
+ M(p(1710),p°) + M(p(1720),p°) + (exchange terms), (33) 

where on the right hand side of Eq. (|33|), the resonances and mesons exchanged in the in- 
termediate states are written inside the bracket explicitly. Each amplitude can be obtained 
straightforwardly by applying the Feynman rules with the relevant interaction Lagrangian den- 
sities. 

For a given invariant collision energy, a/s, the total cross section, a(pp — > pAK + ), can be 
calculated by 

a(pp - pAK+) = jj\M\ 2 *\vr +p 2 - Pa - PA -p K )^t^ (34) 
with the flux factor 

F = 2X 1 / 2 (s,m 2 p ,m 2 p )(2n) 5 , 
\(x,y,z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx, 

where the subscripts 1 and 2 stand for the initial protons, and 3 for the final proton, corre- 
sponding to the diagrams in Fig. [I] for a(BiB 2 — > B 3 YK). This notation, corresponding to 



7 



Fig. |T|, will be used hereafter. \Ai\ 2 is the square of the scattering amplitude, averaged over the 
initial spins and summed over the final spins. To perform the integration over the final state 
three-body phase space, we first integrate on d 3 p 3 and d 3 p\ in the center-of-mass (cm.) frame 
of the final state proton and A, and then integrate on d 3 px in the cm. frame of the initial 
protons |)5] using the Lorentz transformation. Note that our treatment neglects all interference 



terms and final state interactions [pi] , M which are important in near-threshold energy region. 



For the effect of the intereference terms in the model, one can find a discussion in Ref. f44}] . 

Fig. ^| illustrates the separate contributions from n, p and 77-exchanges to the total cross 
sections for the pp — > pAK + reaction, y/s is the pp invariant collision energy in their cm. 
system, while ^/so = m N + m\ + m K is the threshold energy with m N , m\ and m K being 
respectively the masses of the nucleon, A and kaon. One can see that pion exchange is dominant 
at energies near the reaction threshold, while p-exchange is dominant at higher energies. 

As was already discussed in Ref. the total cross section is rather sensitive to the 
value of the cut-off parameter, A p , in the pNN vertex form factor at higher energies. The 
dependence on the values of the cut-off parameter, A p , and coupling constant, g P NN, is shown 
in Fig. |[ After fixing the cut-off parameter to a specific value, the sensitivity of the total cross 
section to the pNN coupling constant, g P NN, (with k = 6.1 for the tensor coupling constant) 
is small. The best values to reproduce the experimental data [49] at energies larger than about 
100 MeV above threshold were found to be, g 2 pNN /^ = 0.74 and A p = 920 MeV. 

Since the total cross sections for the reaction at energies just above threshold receive their 
dominant contribution from pion exchange, we will test the sensitivity of the results to the value 
of the cut-off parameter in the wNN form factor, A^at- Although the cut-off and coupling 
constant, g n NN and A n NN, are not independent in a strict sense, when they are determined 
from the AW phase shift data, we will show the results calculated using two different values 
for the cut-off parameter, A wNN = 3000 MeV and 1050 MeV, with the fixed, usually accepted 
value, g^. NN /4ir = 14.4. The adopted value from this reaction for the cut-off parameter is, 
A wNN = 1050 MeV. 

In Fig. [| we show the total cross sections for the pp — > pAK + reaction calculated with 
the two different cut-offs, A^^n = 3000 MeV (the solid line) and A^attv = 1050 MeV (the 
dashed line), together with the experimental data [[34], |49|]. Results are rather insensitive to the 
value, A n NN- The new data point from COSY |!4j does not seem to be achieved by varying the 
value of the A^nn within a reasonable range. However, it is worth noting that the parameter 
set fixed by the pp — > pAK + reaction, including the cut-off A^nn, must be used also for the 
calculation of the AW — > NYiK reactions. It is not trivial at all if the same parameter set 
could also reproduce the experimental data for the NN — > NUK reactions. Furthermore, 
although the parameters were optimised so as to reproduce the data points around 1 GeV 
above threshold, the recent data for the total cross sections measured at energies about 50 and 
150 MeV above the threshold |47J show an excellent agreement with the predicted values of 



the resonance model [fl4j] . This provides some confirmation of the validity of the parameters 
determined in the model. 

Next, we will discuss the pp — > NY,K reactions. In Fig. |^ we show the separate contri- 
butions from 7r, p and 77-exchange to the total cross sections for the pp — > pH°K + reaction. 
The p-exchange is again dominant at higher energies, while the rj and 7r-exchanges are dom- 
inant near the threshold. The larger contribution from the //-exchange compared to that of 
7r-exchange, which contrasts with the pp — > pAK + reaction, is due to the large ?7AW(1710) 
coupling constant, i.e., the large branching ratio of the A^(1710) resonance to the rjN channel. 
(See also Tables ^ and |2|.) 
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The energy dependence of the total cross sections for the pp — > pT* + K°, pp — > pTPK + 
and pp — > reactions is shown in Fig. |6], together with the experimental data |4"9"| . The 

threshold is, = witv + m s + m ftr ; with the mass of the E-hyperon. The solid and 
dashed lines indicate the results calculated using the two different values the cut-off parameter, 
A-nNN — 3000 MeV and 1050 MeV, respectively. Although the deviations of the experimental 
data are relatively large, our model with the adopted cut-off parameter value, A wNN = 1050 
MeV, reproduces the data fairly well using the values, g 2 pNN /ATr = 0.74 and A p = 920 MeV, 
fixed in the pp — > pAK + reaction. 

In Fig. 0, we show energy dependence of the total cross sections for the pn —* pAK°, 
np — > pTfK and np — > pH~K + reactions together with the experimental data [ 49]j . The 



results are again shown for the two values of A^a^ where the results calculated with the value, 
A n NN = 1050 MeV (the dashed lines in Fig. 0) should be compared with the data. At energies 
about 1 GeV above the threshold, the magnitudes of the total cross sections calculated using 
the two different values of the cut-off, A^nn, become almost equal. This is because, at these 
energies, the momentum of the exchanged pion becomes large and the 7rNN form factor p3 is 
insensitive to the value of the cut-off parameter, A nNN . The model results (the dashed lines) 
overestimate the experimental data roughly a factor 2, or at most a factor of 7 for one data point, 
although the same parameter set reproduces the data for the pp — > pAK + reaction fairly well. 
However, the overestimate of about a factor 2 for the pn — > pAK° reaction, can be contrasted 



with the finding of Faldt and Wilkin | 5T| , that the results for the reaction with a neutron target 
differs by factors of 5 to 10. Their calculation included the iV(1650) resonance alone with pion 
exchange, which is included as one of the contributions in the present calculation. 

We should comment about the discrepancies between the calculated results and the neu- 
tron data. Our calculations for the pn(np) — > NYK reactions are based on SU(2) symmetry 
which should be very good for the present purpose. For the given Lagrangian densities, which 
are explained in Sec. 2, the calculations for the pn{np) — > NYK reactions are trivially related 
to the pp — > NYK reactions through isospin symmetry. We do not have any new mechanism 
to introduce in the model in calculating the pn(np) — > NYK reactions - except for charge 
symmetry breaking, which cannot be expected to lead to such a huge difference between the 
pn(np) — > NYK and pp — > NYK reactions. We should also comment that the data for the 



pn(np) — > NYK reactions tabulated in Ref. JI| are taken from Ref. |50[, whose data contain 



a systematic uncertainty of ~ 13 % in the neutron beam normalization, where that error is 
not included in Fig. [7]. Furthermore, their analysis |S0| was based on a number of hypotheses, 



which in principle, could be strongly model dependent. Thus, we would like to emphasize that 

it is necessary to investigate further the reactions involving a neutron in the initial state, both 
theoretically and experimentally. 

The total cross sections for the K° production channels are obtained by 

cr(pp — > pAK + ) = a{nn — > nAK°), (35) 

a(pn -»• nAK + ) = a(np -> pAK°), (36) 

a{pp -+pY?K + ) = a(nn -> nTPK ), (37) 

a(pp -> nZ+K + ) = a(nn^pZ-K°), (38) 

a{pn -► nY?K + ) = a(np -* pS° <K°), (39) 

a(np -> p^K + ) = a(pn -> nS + K°), (40) 

a{nn -> n^K + ) = a{pp -> pS+K ). (41) 
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3.2 NN AY# 



One expects that contributions from the A" A" — > AYK reactions to the kaon yield in heavy ion 
collisions are small, because of the high threshold energy. Thus, usually these reactions are not 
included in simulation codes. However, there has been no explicit theoretical estimate for the 
reactions. 

In Fig. |8| we show the energy dependence of the total cross sections for the nn — > A~AK + 
and pp — > A ++ £~ K + reactions. At NN invariant collision energies about 3.1 GeV, the mag- 
nitudes of the total cross sections for both reactions become about 70/ifc. These magnitudes 
are comparable with that of the pp — ► pAK + reaction at these energies. However, they become 
about 10 times larger than those of the pp — * pAK + reaction at A^ A" invariant collision energies 
about 3.5 GeV. Thus, these reactions, which are usually discarded in theoretical studies of kaon 
production in heavy ion collisions, might give large contributions to the kaon yield, if there are 
plenty of A^ A" pairs which have such collision energies. 

In order to make an estimate of the uncertainty in the calculation for these reactions, 
we make a comparison with the data for the pp — > NttYK reactions in Fig. || Note that 
these reactions have different thresholds, but the comparison should be helpful in judging 
whether the results contain more than a order of magnitude uncertainties. We restrict ourselves 
to NN invariant collision energies below about 3.5 GeV (as we do throughout this paper) 
the predictions are certainly of the same order of magnitude as the data and therefore not 
unreasonable. We conclude that the NN — > AYK reactions may be equally as important 
as the A^A" — > NYK reactions for heavy ion simulations. However, this argument should be 
checked eventually by heavy ion simulations, keeping in mind the theoretical uncertainties. 

The total cross sections for the other isospin channels are obtained by 



a[nn - 
= 3a (pn 
= 3a (pp - 

a(pp — ► 
: 2a(nn - 
3a (np - 
3a(nn - 
- 6a (pp - 
= 6a (pn 



A~AK+ 
. A°AK + 
A + AK + 

A ++ Z-K + 
-> A-Z°K + 
+ A + Y>-K + 
-> A°E-K + 
-> A + Z°K + 
-> A°E°X+ 



a(pp — 
3a (np 
3a (nn 

a(nn - 
2a (pp - 
3a (pp - 
3a (pn 
6a (np 
6a(nn 



A ++ AK°) 
-> A + AK°) 
-> A°AK°), 

■ a-j: + k°) 

-> A+S+K ) 
-> A°Z + K°) 
-> A + S°K°) 
-> A°S°K°). 



(42) 



(43) 



3.3 AN NYK 



According to the simulation results 
mostly produced through the AA^ 



10, 39, 40, 41, 42 



it is widely believed that kaons are 
NYK reactions in heavy ion collisions. 



NYK and AA 

However, in most of the theoretical studies with the Boltzmann/Vlasov-Uehling-Uhlenbeck ap 
proach (BUU/VUU) 0, [53|, |54], ^6|, 0, or quantum molecular dynamics (QMD) 



39 



the parametrizations used for these total cross sections have been made in a simple phenomeno- 
logical manner in free space || |26], |27 |. In addition, the parametrizations for all isospin chan- 
nels are obtained by fitting to the existing limited isospin channels of the data. Further- 
more, those parametrizations for the reactions involving a A or A 's used to draw the conclu- 
sion |[0, EM HO, ffTl R3, KDf , are obtained using the cross sections for the nucleon except for the 
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isospin difference ||26|| . In addition, the coupling constants relevant for the A are assumed to 
be equal to those for the nucleon, e.g., F NNn = F N ^ W and F^ n = F NNn . 

Recently, this isospin symmetry ansatz suggested in Ref. |26| has been reconsidered Li 
and Ko ||42|| , and elementary kaon production cross sections in baryon baryon reactions have 
been calculated using a one pion and one kaon exchange model. But their approach is still 
not fully consistent because they adopted the results of the pion baryon reactions from the 



resonance model |4q] , where the resonance model [[lg] does not include one kaon exchange in 
its basic assumptions. Thus, if one wants to be consistent with the parameters determined in 
the pion baryon reactions, one kaon exchange should not be introduced. 

In Fig. we show the energy dependence of the total cross sections for the A ++ n — > 
pAK + and A~p — > n12~K + reactions, together with the data for the pp — > pAK + , pp — > pJ2 + K° 
and pp — > riL + K + reactions |49| plotted at the same excess energy above each threshold energy 
corresponding to the A and E production channels. The magnitude of the total cross section 
for the A ++ n — > pAK + reaction is almost equal to that of the pp — > pAK + reaction at the same 
excess energy just above the threshold, \/s—*/sa < 100 MeV. However, at energies ^/s — ^/so ~ 1 
GeV, the magnitudes of both the A 



++ 



n 



pAK + and A~p — > nYj~K + reactions become larger 
by about a factor of 6 than those of the pp — > pAK + , and pp — > pH + K° and pp — > n'E + K + 
reactions, respectively. One of the main reasons for this is the different spin structure of nucleon 
and A, as can be seen from Eq. (|3l|) , which produces different energy dependence in the Lorentz 
invariant scattering amplitude. (The difference in form factors and coupling constants is also 
the other reason.) Our result may be consistent with the conclusion that the AiV — > NYK 
reactions give the largest contribution to kaon yield in heavy ion collisions. However, this 
should be explicitly checked. 

The total cross sections for the other isospin channels of the AiV — > NYK reactions are 
obtained by: 



a(A ++ n - 


-> pAK + ) 


= o(A~p — 


> nAK°) 


3a(A + j9 - 


-> pAK + ) 


= 3a(A + n - 


-> pAK°) 


3<x(A+n - 


-> nAK + ) 


= 3a(A°p- 


-> pAK°) 


3a(A°p - 


-> nAK + ) 


= 3a{A°n- 


-> nAK°), 



a(A~p - 
2a(A ++ n 
= 3cr(A + n - 
= 3a(A°p - 
= 3a(A°n - 
= 6a(A + p 
= Qa(A + n 
= 6<r(A°p 



nT,-K + 

> pY?K + 
pT,-K + 
pY>-K + 

> pY?K + 

> nY?K + 

> nY?K + 



n 



a(A+ 
2cr(A~p 
3a(A + p 
3a(A + n 
3a(A°p 
Qa(A + n 
6a(A°p 
6(x(A n 



+ p j: + k°) 

> nY?K G ) 

* pZ°K ) 
pY?K°) 



(44) 



(45) 



3.4 AiV -> AYK 

Similarly to the situation for the NN — > AYK reactions, the AA — > AYK reactions are not 
usually included in the simulation codes due to the high threshold energy. However, these 



reactions may be responsible for the low energy part of the kaon spectra. In Figs. |TT] - [14], we 



11 



show energy dependence of the total cross sections sufficient for obtaining all isospin channels 
in the reactions, together with the experimental data of the pp — > pAK + , pp — > pH°K + , 
pp — > p£ + if° and — > reactions plotted at the same excess energies above the 

corresponding thresholds. The magnitudes of the total cross sections are almost equal to those 
measured for the pp — > NYK reactions at the same excess energy, — yfso > 100 MeV, with 
respect to each corresponding threshold, ^/sq. 

The total cross sections for the other isospin channels are obtained by: 

a(A ++ p -> A++A^ + ) = a(A ++ n — > A ++ AK°) 

= a(A~p -> A~A^ + ) = o-(A _ n -> A"AK°), (46) 



3(7(A + n -> A°A^+) = 3a(A°p -> A+A^°) 

= 4a(A ++ n -> A+A^+) = 4o-(A + p — > A ++ AK°) 

= 4a(A°n A-A7T+) = 4<r(A-p -> A Aif°), (47) 

(j(A + p -> A+A^+) = a(A + n -> A+A^°) 

= o-(A°p — > A°Alf + ) = <r(A°n -> A°AiT ), (48) 

( j(A ++ n -> A ++ S-^+) = <r(A ++ p -> A++S+K ) 

= cr(A-n -> A-E-X+) = a{A-p^ A-£ + K°), (49) 

3cr(A p -> A+£-^+) = 3 ( x(A + n -> A°S + K°) 

= 4a(A + p -> A++S-^+) = 4 ( j(A ++ n -> A+S+K ) 

= 4cr(A~p -> A°S-^+) = 4 ( j(A°n ^ A-S+^°), (50) 

^(A+n -> A+S-^+) = ^(A+p -> A+£+^°) 

= a(A°n -> A°£-^ + ) = a(A°p -> A°S+^°), (51) 

<t(A ++ p -> A ++ S°^+) = a(A ++ n -> A ++ S°^°) 

= (j(A-p -> A-S°K + ) = tr(A-n -> A-£°K°), (52) 

3<7(A+n -> A°S°^+) = 3<r(A°p -> A + S°K°) 

= Aa(A ++ n -> A + S°K + ) = 4a(A + p -> A ++ £°^°) 

= 4a(A°n ^ A-S°^+) = 4a(A-p -> A°S°K°), (53) 

^(A+p -> A + E°^ + ) = ^(A+n -> A + S°K°) 

= a(A°p -> A°£°^ + ) = ^(A ™ -> A°S°K°), (54) 

3a(A + p -> A°£ + ^ + ) = 3a(A°n -> A+S^K ) 

= 4a(A++p -> A+S+ir + ) = 4a(A+n ^ A++S-K ) 

= 4cr(A°p -> A-Z+K+) = 4a(A-p^ A°Z°K°). (55) 
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3.5 AA AYK 



The energy dependence of the total cross sections for the AA 
reactions are shown in Figs. 15, [16] and |T 



AAK and AA -> AEif 
together with the experimental data for the pp — > 
pAK + and — > pTPK + reactions plotted at the same excess energies above the corresponding 
hyperon production threshold. As expected, the magnitudes of the total cross sections for these 
reactions are small, because of the high threshold energies and three A's are involved in the 
processes. The magnitudes for these reactions are almost 10 times smaller than those for the 
pp — > pAK + reaction which is shown in Fig. |17. The contributions from these reactions to the 



kaon yield in heavy ion collisions are therefore negligible. 

The total cross sections for the other isospin channels are obtained by: 



a(A+A++ - 


-> A ++ AK + ) 




a(A°A- - 


A~AK°) 


= 3<r(A°A° 


-> A~AK + ) 


— 


3a(A+A+ 


-> A ++ AK°) 


= 9a(A+A+ 


-> A+AK+) 





9a(A°A° - 


■> A°Aif°), 


a(A°A++ 


-> A + AK + ) 





a(A + A- - 


> A°AK°), 


ct(A°A + 


-> A°Afs: + ) 




a(A+A° - 


A + AfT°), 


a (A++A° -> 


A ++ E-^ + ) 


= 


a(A ++ A° - 


-> A ++ S°K°) 


= ct(A ++ A~ - 


■>■ A + E~^ + ) 




a(A ++ A- 


-> A+S ^ ) 


= a(A-A ++ 


-> A £°tf+) 




a(A-A ++ 


-> A S + iT°) 


= a(A~A + - 


-> A-S°K + ) 


= 


a(A~A+ - 


♦ A-£ + K°), 


a(A-A° - 


-> A-X-K+) 




a(A++A+ 


-> A ++ E+^° 


= 2a(A + A ++ -» 


A ++Z°K + ) 




2a{A°A- - 


-> A-S°K°) 


= 3a(A + A + -> 


A ++ E~K + ) 




3a(A°A° - 


-> A-E+iT ) 


= 6<r(A°A° - 


-> A-S°iT + ) 




6a(A+A+ 


-> A ++ S°K°) 


= 9a(A°A° - 


-> A S-^ + ) 




9a(A + A+ 


-> A+S+K ) 


= 18<r(A + A+ - 


-> A + E°K+) 




18a(A°A° 


-> A°S°K ), 


2a(A°A++ - 


-> A + E°K + ) 




2ct(A + A~ 


-> A°S°K°) 


= 3<r(A°A + - 


■> A+ITK+) 




3a(A+A° - 


-> A°E + K°), 


a(A-A+ - 


-> A £-AT + ) 




a(A++A° - 


-> A+E+if ) 


= 6<r(A°A + 


-> A°S°K + ) 




6a(A+A° - 


-> A+E°K°). 



(56) 
(57) 
(58) 



(59) 



(60) 



(61) 



(62) 
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4 Parametrizations 



The energy dependence of the total cross sections calculated so far will be parametrized in the 
following form: 

a(B 1 B 2 ->B 3 YK) = a(j--l) ^J, (63) 

where s and sq are respectively the squares of the invariant collision energy and the threshold 
energy, a, b and c are the parameters to be determined so as to reproduce the calculated energy 
dependence of the total cross sections, with a in units of mb. 

As was demonstrated in Ref. plj , the functional form, Eq. (|63|), can reproduce quite 



well the near-threshold behaviour of the total cross sections. The reason is that the function 
reflects the energy dependence of the phase space reasonably just above the threshold, where 
the energy dependence of the amplitude is usually weak. In addition the particles in the final 
state may be treated nonrelativistically due to the small amount of energy available. Assuming 
the relevant amplitude squared, |T| 2 , to be a constant, the total cross section can be written as 

o = f |T| 2 , (64) 

with i? 3 the three-body phase space volume and F the flux factor. In Fig. [18| we show the results 
for the energy dependence of the total cross section for the pp — > pAK + reaction, calculated 



using Eq. (p4|), together with the experimental data |3J, f49| . The square of the constant 
amplitude, |T| 2 , was determined so as to fit the data point of COSY |3"4]] . The parameter, b, 
which is mostly relevant to reproduce the shape of the energy dependence of the total cross 
section near the threshold, is obtained to be b = 1.995 in this case. This value of b gives a 
behaviour of the energy dependence very close to that of phase space in the nonrelativistic 
|66]| . On the other hand, the parameter, b, obtained to reproduce our results is, b ~ 2.000, 



case 

which is very close to the value of the above phase space calculation. The last term in Eq. 
(]63]) ) (s /s) c , reflects both the energy dependence of the flux factor which becomes important 
at energies y/s — ^/sq > 0.5 GeV, and also the contribution from meson exchange which is 
expected to behave as s~ 2 . 

The parameters, a, b and c in Eq. fl5"3"D, determined so as to reproduce the results calculated 
in our model, are listed in Table [| They are sufficient to reproduce all the isospin channels 
treated in the previous section. The corresponding threshold energies squared, so, are also 
indicated. 

Note that the parameters given in Table § are determined so as to reproduce the theo- 
retical results at energies up to about 2 GeV above the corresponding threshold energies. The 
reason is, when kaon production is studied at energies larger than 2 GeV above threshold, we 
need to include more resonances which decay to kaon and hyperon in the model, as well as tak- 
ing into account the other kaon production processes with more than three particles appearing 
in the final states, where it seems impossible due to the present experimental data available. 
In practice, we believe that the parametrizations are recommended to be used up to invariant 
collision energies about 3.6 GeV, which are usually enough. 

Although we have parametrized explicitly all the isospin channels required, it has been 
traditional to use isospin-averaged cross sections. In order to compare our results with those 
isospin-averaged cross sections, we have also calculated them within the model. The results 
are again parametrized by the function, Eq. (^). The parameters, a,b and c, for the isospin- 
averaged cross sections, are listed in Table f|. Note that the parameterizations in Table |4] also 
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Table 3: Parameterizations for energy dependence of the total cross sections, a, b and c, 
appearing in Eq. ( |63"D were determined so as to reproduce the calculated total cross sections 
at energies, yfs — J~Sq < 2 GeV, with yfs and ^/so being, respectively, the invariant collision 
energy and threshold energy. The parametrizations are recommended to be used up to invariant 
collision energy about 3.6 GeV. 



No. 


Reaction 


so (GeV 2 ) 


a (mb) 


b 


c 


1 


pp — > pAK + 


6.504 


1.879 


2.176 


5.264 


2 


pn — ► nAK + 


6.504 


2.812 


2.121 


4.893 


3 


pp — > pH°K + 


6.904 


5.321 


2.753 


8.510 


4 


nn — > n£~.K'+ 


6.904 


7.079 


2.760 


8.164 


5 


— > riE°K + 


6.904 


6.310 


2.773 


7.820 


6 


np — > pY,~K + 


6.904 


11.02 


2.782 


7.674 


7 


pp — ■> 


6.904 


1.466 


2.743 


3.271 


8 


nn -> A _ A^+ 


8.085 


6.166 


2.842 


1.960 


9 




8.531 


10.00 


2.874 


2.543 


10 


A++n -> pAfT+ 


6.504 


8.337 


2.227 


2.511 


11 


A~p -> n£-ir+ 


6.904 


52.72 


2.799 


6.303 


12 


» J |_ » _| 1 i T r 1 

A++p -> A++A-FT+ 


8.085 


2.704 


2.303 


5.551 


13 


A+n -> A°AK+ 


8.085 


0.312 


2.110 


2.165 


14 


A+p -> A+A/ST+ 


8.085 


2.917 


2.350 


6.557 


15 


A++n -> A++E-AT+ 


8.531 


10.33 


2.743 


8.915 


16 


A°p A+S-iT+ 


8.531 


2.128 


2.843 


5.986 


17 


A+n -> A+S-ir + 


8.531 


10.57 


2.757 


10.11 


18 


A++p -> A++£°K+ 


8.531 


10.30 


2.748 


9.321 


19 


A+n -> A°S°ir+ 


8.531 


1.112 


2.846 


5.943 


20 


A+p -> A+£°if+ 


8.531 


10.62 


2.759 


10.20 


21 


A+p -> A°S+K+ 


8.531 


0.647 


2.830 


3.862 


22 


A+A++ -> A++AK+ 


8.085 


1.054 


2.149 


7.969 


23 


A°A++ -> A+AK+ 


8.085 


0.881 


2.150 


7.977 


24 


A°A+ -> A°AK+ 


8.085 


0.291 


2.148 


7.934 


25 


A++A -> A++S-iT+ 


8.531 


3.532 


2.953 


12.06 


26 


A-A° -> A-S-iT+ 


8.531 


7.047 


2.952 


12.05 


27 


A°A++ -> A+£°^+ 


8.531 


2.931 


2.952 


12.03 


28 


A"A+ -> A°E-i<:+ 


8.531 


5.861 


2.952 


12.04 
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Table 4: Parameterizations for the energy dependence of the isospin-averaged total cross 
sections. (See also the caption of Table 0.) For the K + production reactions, replace the values 
for the parameter a by ha. 



No. 



Reaction 



a (mb) 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 



NN 
NN 
NN 
NN 
AN 
AN 
AN 
AN 
AA 
AA ■ 



■ NAK 

■ NZK 
• AAK 

AEK 
NAK 
NTjK 
AAK 
AT,K 
AAK 
AY.K 



2.330 
15.49 
9.249 
12.50 
4.169 
39.54 
2.679 
21.18 
0.337 
5.140 



2.140 
2.768 
2.842 
2.874 
2.227 
2.799 
2.280 
2.743 
2.149 
2.952 



5.024 
7.222 
1.960 
2.543 
2.511 
6.303 
5.086 
8.407 
7.967 
12.05 



include K° production channels. Thus, the parametrizations for the K + production channels 
alone can be obtained replacing the values for the parameter a in Table f| by \a. 

We show the energy dependence of the isospin-averaged total cross sections, a, calculated 



in the model for all the baryon baryon reactions in Fig. |T^. Note that the horizontal axes indicate 
invariant collision energies. Our results illustrate that the magnitudes of the isospin-averaged 
cross sections at energies, ^fs - < 200 MeV, are, a(AN -> NY K) ~ a(NN -> NYK) 
for both the A and X production channels. However, at energies, \/s — t/sq > 200 MeV, the 
magnitudes of the cross sections a(AN — > NYK) become about 6 times larger than those 
of a(NN —>■ AYK). It is also noticable that the magnitudes of a(NN — > AYK) become 
comparable with those of a(NN —>■ NYK) at NN invariant energies about 3.1 GeV, and 
overcome factors of 5 to 10 at iViV invariant energies about 3.5 GeV. In particular, this is 
pronounced for the A production reactions. It is not clear at the moment how much the 
NN — > AYK reactions contribute to the total kaon yield in heavy ion collisions because 
they have usually been neglected. It seems worthwhile to perform calculations including these 
reactions. At invariant energies larger than about 3.2 GeV, the magnitudes of a(AN — > AYK) 
are also comparable to those of a(NN — > NYK). 

It is interesting to see if the relation for the elementary kaon production cross sections, 
a(NN -> AAK) > a(AN -> NAK) > a(NN -> NAK), can appreciable influence inclusive 
kaon production in proton-nucleus or heavy ion collisions. If this is the case, high energy 
inclusive kaon production in heavy ion collisions will be also significantly influenced by the A. 
A test for this argument can be made by performing calculations for high energy inclusive kaon 
production both, neglecting all the A's, and including them, and comparing the results with 
the existing experimental data. 

Next, we compare the isospin-averaged total cross sections calculated in this model, with 
the commonly used relation suggested by Randrup and Ko p6f . For the NN — > NYK reactions, 
they suggested the relations: 

a(NN -> NAK) = 3 a(pp -> pAK + ) } (65) 
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a(NN -> NEK 



a(pp -> pS°ir + ) + <r(pp -> pE+AT ) 



(66) 



In Fig. we show the isospin-averaged total cross sections calculated in the model, 
a(NN — > NAK) and a(NN — > NEK) (denoted by our model), and the quantities calculated 
using the right hand side of Eqs. flBDp and ( pB| ) (denoted by prescription from [26]). The explicit 
calculation shows disagreement of about a factor 2 between the relation suggested by Randrup 
and Ko [26] for most of the energy range up to 2 GeV above the thresholds, for both the A 
and E production reactions. Thus, the calculations performed using the relation suggested 
by Randrup and Ko |2f| to obtain the isospin-averaged cross sections for the iViV — > NYK 
reactions, may have about a factor of 2 ambiguity for these reactions. 

As for the A A" — > NYK reactions, they suggested the relations |^6| : 



a(AA-> NAK) 
a(AN -> NEK) 



9 

- a(pp -> pAK + ), 
9 



4 



a(pp -»■ pS AT + ) + a(pp -> pE + K°) 



(67) 
(68) 



In Fig. we show similar quantities calculated within the present model as for the 
A^ A" — > NYK reactions. At energies near the thresholds, y/s — < 200 MeV, our results 
are in good agreement with the relations Eqs. ( |67"D and fl68"D , suggested to obtain the isospin- 
averaged cross sections. However, at energies, yfs — ^/so ~ 1 GeV, our results show a factor of 
6 discrepancy with the relations suggested by Randrup and Ko |26| . 

It is worth noting that the recent calculation for kaon production in heavy ion collisions 
in free space |[35|| , performed at energies near or below the threshold region of the A^A" — > NAK 
reaction, illustrates that the dominant contribution comes from the secondary, pion induced 
Cassing et al. [ 



reactions. Cassing et al. [|3J| and Li and Ko |]42] used the relations of Eqs. (|65| ) - (|68|), to 
obtain the isospin-averaged cross sections. In view of our results, it seems necessary to perform 
calculations using the total cross sections parametrized in a consistent manner for both the 
baryon baryon and pion baryon reactions. 



5 Application to heavy-ion collisions 

Kaon production in heavy- ion collisions at SIS energies (1-2 AGeV) has been studied theoreti- 
cally using various transport models. Important inputs for such studies are the elementary kaon 
production cross sections in pion baryon and baryon baryon interactions. Conclusions concern- 
ing nuclear matter equations of state and kaon in-medium properties, obtained by comparing 
calculations with measurements, depend on the input kaon production cross sections. 

In Ref. |n|, a detailed study of these elementary cross sections was carried out, based 
on both the parametrization of experimental data and theoretical model calculations. Here we 
want to compare the kaon yields in heavy-ion collisions obtained using the cross sections in this 
work with the cross sections of Ref. ]19| . We will use the same relativistic transport model as 
in Ref. [TP| so that the difference in kaon yield is due entirely to the different elementary cross 



sections used. In this way, we can get a feeling for the sensitivity of transport model predictions 
to the underlying elementary cross sections. 

We take central Au+Au collsions at 1 AGeV and Ni+Ni collisions at 1.8 AGeV as exam- 



ples. The results are shown in Figs. |22] and respectively. From these figures, we can make 
the following observations: 
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Generally, the kaon yield obtained with the cross sections from this work is smaller than 



that obtained with cross sections from Ref. ||19|| . The difference is about a factor of two 
at 1 AGeV, and reduces to about 40% at 1.8 AGeV. 

Using both sets of cross sections, the A A channel is found to make the largest contribu- 
tions among all the baryon baryon interactions. This is also in agreement with the earlier 
findings using the simple Randrup-Ko parametrization. 

The contribution from the AA channel is particularly small using the cross sections 
from this work. This is due to the fact that in the present model, the cross section for 



A A — > NYK channel is zero, while in Ref. (T^j off-shell pion exchange was included and 
gave a small contribution. 



6 Summary 

Within the resonance model, which is fully relativistic, but includes only the tree-level processes 
as an effective and empirical contributions, we have studied the elementary kaon production 
reactions in baryon baryon collisions in a consistent manner with those for the pion baryon 
reactions. We have calculated and parametrized the energy dependence of the total cross 
sections for all the necessary isospin channels in the AA — > NYK, AA — > AYK, AN — > 
NYK, AN -> AYK and A A -> AYK reactions (Y = A and E). We have also calculated and 
parametrized those for the isospin-averaged cross sections and compared with the the relation 
suggested by Rundrup and Ko to obtain the isospin-averaged cross sections. 

The present study has several points which may need to be improved, for example, con- 
sistency with the unitarity, the possibility to add other types of meson exchanges, such as kaon 
exchange, and the overestimate for the pn(np) — > NYK reactions and the reactions involving 
the A, we would like to emphasize a practical aspect of this work, which we belive providing 
a useful inputs for further investigations of kaon production in heavy ion collisions. Keeping 
these in mind, we state a few comments as follows. First, our explicit calculation indicates 
that the contributions from the AA — > NYK and AA — > AYN reactions are substantially 
larger than those from the AA — > NYK reactions at higher energies. Second, the relation 
suggested by Randrup and Ko |26|] (c.f. Eqs. ( pop and (pop) deviates from our results by about 
a factor of 2 for the AA — > NYK reactions. For the AA — > NYK reactions, the relation (c.f. 
Eqs. ( |o7| ) and (p8"P) holds well at low energies, up to about 200 MeV above the threshold, but 
deviations by as much as a factor of 6 have emerged at higher energies. From these facts we 
conclude that it is necessary to use consistent parametrizations of total kaon production cross 
sections for both the pion baryon and baryon baryon reactions. Furthermore, one should use 
isospin-averaged total cross sections calculated in a consistent manner. Lastly although the 
present study has been made for exclusive kaon production, there might also be a significant 
difference between the AA and AA collision reactions in inclusive kaon production in proton- 
nucleus or heavy ion collisions. Thus, it is interesting to perform calculations for inclusive kaon 
production, focusing on the effect of the A. 

Acknowledgement 

We thank J. Haidenbauer for helpful comments on the unitarity and coupled channel calcula- 
tion. K.T. would like to thank K. Yazaki for valuable discussions at the very early stage of this 
work. A.S. would like to thank W. Cassing and U. Mosel for valuable discussions. This work 



18 



was supported in part by the Forschungszentrum Jiilich, the Australian Research Council, and 
the U.S. Department of Energy under grant No. DE-FG02-88ER40388. 



7 Appendix 

7.1 Coupling constants of the resonances 

Here we give relations between the branching ratios of the adopted resonances and the corre- 
sponding coupling constants squared, evaluated in the rest frame of each resonance. Note that 
all the coupling constants appearing below should be understood as being multiplied by the 
appropriate form factors. In addition, the confidence level, spin and parity of each resonance 
is specified. For the definition of \{x, y, z) appearing below, see the text, after Eq. 



7.1.1 iV(1650)****(i; 



r(iV(1650) -> Ntt) 

with \p N \ 
r(iV(1650) -> Nrj) 

with \p N \ 
r(iV(1650) -> Att) 

with |pa I 
r(iV(1650) -> AK) 

with |pa| 



7.1.2 N(1710)***U) + 



^Inn (1650) {En + m N ] 
4vr 77Jat(1650) 



\P 



Xl/ \ m N(1650)i m Ni m l< 



2m 



Af(1650) 



(69) 



9% NN (1650) {En + win] 



Att 



m iV(1650) 



■I/ 7 



A 1/2 (m^ (1650) ,m^,m^) 



2m 



iV(1650) 



, 9%AN(1650) fnN(1650){El 



Qtt 



* 1/2 W(i650)» m A»ra£) 

2mAT(i 650 ) 
9\kN(1650) {E\ + 771a) 
7™jV(1650) 



47T 



PA\ 



A 1/2 (m^ (1650) ,mj,m^) 

2mAr(i 650 ) 



(70) 



\Vt 



(71) 



(72) 



T(iV(1710) -> A%) 
with |/j/v| 



o 9ttNN(i710) {En — tun) , _ , 

d ~ A \Pn\, 

4tt mjv(i7io) 

* 1/2 (™N(1710), m N,™l) 

2mAT(i7i ) 



(73) 
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r(7V(1710) -> Nrj) 

with |ptv I 
r(7V(1710) -> Np) 

with \pn\ 
r(JV(1710) -> Att) 

with |pa| 
r(7V(1710) -> AX) 

with |pa| 

r(Jv(i7io) -> sx) 

with 

.1.3 iV(1720)****(§) + 

r(iV(1720) -> A%) 
with |pjv| 
r(iV(1720) -> iVr/) 

with |pjv| 
r(iV(1720) -> Np) 

with |pjv| 
r(JV(1720) -> Att) 

with |pa. I 



9%nn(i7io) {En — it^n) 



^JV(1710) 



\Pn\, 



Al/2 K(mo)- m M. m ?) 



2m 



#(1710) 



? gp##(1710) |p 



'# 



47T 



"T#(1710) 



mf 



A 1/2 i 



m #(1710)> m #> m p) 



2mAr(i7i ) 

O fi f2 -A#(1710) ^JV(1710) (Ea + ^a) i -. 



677 



2 



■Pa 



m #(1710)> m L 



2m 



#(1710) 



gLr#(mo) (£a - m A ) _ 

47T m A r (1710 ) 

Al/2 ( m Af(mo)^ m A^ m i-) 



2m 



= 3 



■#(1710) 

fi , SA'Af(1710) (■E's — ™e) 
477 mAT(i7i ) 

Al/2 (^AT(1710)> TO L^) 



IpsI, 



2m 



#(1710) 



'AT 



g ff 2 ##(1720) (£# + mjy) ^ 
1277 m N (i 7 2o)ml 

Xl/ H m N(1720)' m N,ml) 

2m A T( 172 o) 

qInn (1720) (En + ra#) . _ 
127T m Ar (i72 )m 2 

A 1/ ' 2 (m 2 V( - 1720 - ) , mfy, m 2 ) 

2m A r( 172 o) 

gflp##(1720) (-E# — ™>n) 
12n m N (172Q) 

Xl/ H m N(1720)' m Ni m2 p) 

2mjv( 1720 ) 

2 f 2 A#(i72o) ^a I Pa I ' ^ E A 
367T mAr(i 72 o) I- ™A 

Xl/2 (^l 2 N(1720), m A^l) 

2m N (i 72 o) 



■N , 



Pn 



3 + ^ 



|p#l, 



m A 
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T(iV(1720) -> AA) 
with 

r(JV(1720) -> £A) 
with |ps| 



7.1.4 A(1920)***(§) + 



9k AN(1720) {E\ + 7T7-A 

127T m N{ i 72 o)m" K 
A 1/2 (m^ (1720) ,m|,m^) 



2 Ipa| 3 , 



2m A r( 172 o) 
„ fl r KKJV(1720) (As + m E ) 3 

127T m N (i 72 o)m K 
X 1/2 (m 2 Nil720) ,ml,m 2 K ) 



2m 



JV(1720) 



T(A(1920) -> JVtt) 
with |pjv| 

T(A(1920) -> £A) 
with |p~s| 



5VjVA(1920) (AjV + ?™iv) | -» 1 3 

12tt m A(192 o)m^ 



2m 



A(1920) 



^7rAfA(1920) (Ag + m S 
127T m A (1920)«% 

AV2 ( m A(1920) ; m L m x) 

2m A (i920) 



2 fel 3 , 



7.2 Cross section relation for the 7riV — > AAT reactions 

Here, we add a relation for the total cross sections for the ttN — * AA reactions, which 
mentioned in Ref. ||46| . The relation is: 



a(n + n 



AA+) 
2a(n°p -> AA+) 



a{n p 



AA C 



2(7(7r n -> A A ). 
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B 3 (p 3 ) Y(p 4 ) K(k) B 3 (p 3 ) Y(p 4 ) K(k) 




B,(p,) B 2 (p 2 ) B,(p,) B 2 (p 2 ) 
a) b) 



Figure 1: Kaon production processes in the resonance model. Bi (i = 1,2,3), Y and R stand 
for respectively, either the nucleon or the A, either the A or the E hyperon, and the baryon 
resonances. 
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Figure 2: Separate contributions from ir, p and 77-exchanges to the total cross section for the 
pp — > pAK + reaction, s 1 / 2 and s l J 2 = itln + rn\ + itlk are the invariant collision energy and 
threshold energy, respectively. 
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Figure 3: Dependence on the pNN coupling constant and cut-off parameter of the total cross 
section for the pp — > pAK + reaction. The dots show experimental data from Ref. [49 . 
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Figure 4: Dependence on the cut-off parameter in the irNN form factor, A^nn, of the total 
cross sections for the pp — > pAK + reaction. The dots and square are the experimental data 



respectively from Ref. [49 and Ref. p 
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Figure 5: Separate contributions from n, p and 77-exchanges to the total cross sections for the 
pp — > pTPK + reaction. 
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A nNN =3000 MeV 



A nNN =1050 MeV 



80 - 




j i i i i_ 



s 1/2 -si /2 (GeV) 



Figure 6: Energy dependence of the total cross sections for the pp — > p~E + K°, pp —>■ pTPK + 
and pp — > n'E + K + reactions. The dots stand for the experimental data from Ref. with 
error bars. The solid and dashed lines show our results calculated using different values for the 
cut-off parameter, A^nn = 3000 MeV and 1050 MeV, respectively. 
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Figure 7: Same as in Fig. || but for the pn —>■ pAK°, np —>■ pTfK and np 
reactions. Error bars are smaler than the dots. 
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s 1/2 -s 1 f (GeV) 

Figure 8: Energy dependence of the total cross sections for the nn — > &~K.K + and 
A ++ S~i ; r + reactions. 
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Figure 9: Comparison of the energy dependence of the total cross sections for the pp — > AYK 
and pp — > NirYK reactions. Note that the parametrizations are recommended to be used up 
to invariant collision energy about 3.6 GeV. 
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Figure 10: Energy dependence of the total cross sections for the A ++ n — > pAK + and A~p 
K + reactions. The dots, triangles and squares are the experimental data for the pp 
pAK + , pp — ► pYi + K° and pp — > nS + K + reactions, respectively, from Ref. {49 . 
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Figure 11: Energy dependence of the total cross sections for the A ++ p — > A ++ AK + , A + n — 
A°AK + and A + p — > A + AK + reactions. The dots are the experimental data from Ref. P5| 
plotted at the same excess energies from each threshold energy. 
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Figure 12: Energy dependence of the total cross sections for the A ++ n — > A ++ £ K + , 
A + T,~K + and A+n -> A+S"^ + reactions. 
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Figure 13: Same as in Fig. |TT], but for the A ++ p -> A ++ £°if + , A+n -> A°S°i^ + and 
A + p — > A + S°i^ + reactions, and the data for the pp — ► pTPK + reaction. 
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Figure 14: Energy dependence of the total cross sections for the A + p — > A°'E + K + reaction 
together with the experimental data from Ref. jy| for the pp — > pH + K° and pp — > nE + K + 
reactions, plotted at the same excess energies. The threshold for the A + p — > A°E + ii' + reaction 
is = m A + ms + m K . 
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Figure 15: Same as in Fig. [11], but for the A + A~ 
A°A+ -> A°AK + reactions. 
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Figure 16: Energy dependence of the total cross sections for the A ++ A° — > A ++ £ K + , 
A- A -> A-^-R + and A~A+ -> A £-fT + reactions. 
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Figure 17: Same as in Fig. [14] but for the A A 
data are for the pp —>■ pTPK + reaction. 



A + 12 K + reaction and the experimental 
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Figure 18: Phase space consideration of the total cross section for the pp — > pAK + reaction, 
according to Eq. (|4]). The dots and square show the experimental data from Ref. f|9| and 
Ref. [Q, respectively. 
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Figure 19: Energy dependence of the isospin-averaged total cross sections for the A and £ 
production reactions in baryon baryon reactions. 
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Averaged cross sections 




Figure 20: Energy dependence of the total cross sections both calculated in the model, isospin- 
averaged, and using the relations Eqs. (|B5D and suggested by Randrup and Ko P^| . The 
solid lines show the model calculations for the isospin-averaged total cross sections (denoted 
by "our result"), while the dashed lines show the model calculations obtained using the right 
hand side of Eqs. (|65| ) and ( j66| ) (denoted by "prescription from p6|"). 
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Figure 21: Same as in Fig. [2(], but for the AiV — > NAK and AiV — > NT,K reactions and 
Eqs. (|67|) and (|68| ), respectively. 
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Figure 22: Kaon yield in Au+Au collisions at 1 AGeV and 1 fm obtained with two sets of 
elementary cross sections. 
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Figure 23: The same as Fig. 22, for Ni+Ni collisions at 1.8 AGeV. 
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